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1 Introduction, SCRMs and Lovasz-Schrijver proce- 
dures 



Since 1960's, complementarity problems attracted a very significant attention in the theory 
as well as applications of operations research. See, for instance, the book on LCP [4]. In this 
paper, we consider various complementarity problems in the context of Successive Convex 
Relaxation Methods (SCRMs) proposed by the authors [5, 6]. Since these methods can be 
used to compute the convex hull of any compact subset of an Euclidean space described by 
a system of quadratic inequalities and a compact convex set, they can be used to attack 
many complementarity problems from several angles. 

In the special case of 0-1 optimization problems over convex sets, or more specially 
polytopes, there are many Successive Convex Relaxation Methods (SCRMs) based on lift- 
and-project techniques. We also discuss some of the relationships of general SCRMs and 
these more specialized algorithms in solving LCPs. 

Let F be a compact set in the n-dimensional Euclidean space R n . SCRMs take as input, 
a compact convex subset Co of R n and a set Vf °f quadratic functions which induce a 
description of F such that 

F = {x E C : qf(x; 7 , q, Q) < 0, qf(-; 7 , q, Q) G V F }. 

Here we denote by gf(-]'y,q,Q), the quadratic function (7 + 2q T x + x T Qx). Note that the 
variable x is irrelevant outside a context and it will always be clear what the variable vector 
is, from the context. 

Let £ be an integer such that 1 < 2£ < m, d G R m , and let A be a compact convex 
subset of R m . Consider the convex optimization problem with complementarity conditions: 



maximize d T u 

subject to u E A, < Ui, < u i+ t, UiU i+ t = 0, Vi G {1, 2, . . . , £}. 



(1) 



First of all, it is clear that LCP, with a known upper bound on a solution of it, is a special 
case of (1) (we can take m = 2£ and A as an affine subspace intersected with a large enough 
ball). Secondly, it is very elementary to formulate this problem as a mixed 0-1 optimization 
problem with convex constraints: 



maximize c v 

subject to v G C , Vi G {0, 1}, G {m + 1, m + 2, . . . , n}, 



(2) 



where 



n 



v 



( u \ 

V m +1 

V v n j 



G R 



m+l 



u G A, 

<Ui< rv m+i , 

< u i+l < r(l - v m+i ), 

Vie {1,2,..., £} 




m + I, r > max {max{t/j : u G A}} . 



G R 



m+l 
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In general, we allow C to be an arbitrary compact convex set in R n . There are various 
successive convex relaxation methods that can be applied to such a problem. 

We can represent the feasible region F C R n of (2) as 

F = {v G Co : p(v) < 0, Vp(-) G Pf}, 
where denotes a set consisting of quadratic functions 

{vj-vi), {-vf + Vi), i G {m + l,m + 2,...,n} 

on R n . 

In connection with the SCRMs and also the Lovasz-Schrijver procedures (see [8]), it 
seems convenient to introduce the following notation: For every compact convex relaxation 
C C C of F and every subset D of D = {d G R n : ||d|| = 1}, 



V 2 (C, D) 
M(C, D) 



N + {C,D) 



{-{d T v - a(C , d))(d T v - a(C, d)):deD, d £ D}, 
3V e S n such that 
v E C : 7 + 2q T v + Q • V < 0, 

Vqf(-n,q,Q) EV F UV 2 (C\D) 

(a Semi-Infinite LP relaxation of F), 

1 v T \ 
x v V ) 
~( + 2q T v + Q*V < 0, 

Vqf(-,l,q,Q)eV F UV 2 (C,D) 

(an SDP relaxation of F), 



3V G <S n such that 



« G Co 



where a(C, d) = max{d T i; : u G C} for every d £ D. Let <S" and <S^ +n denote the set of 
n x n symmetric matrices and the set of (1 + n) x (1 + n) symmetric positive semidefinite 
matrices, respectively. The corresponding variants of Successive Semi-Infinite LP Relaxation 
Method (SSILPRM) and Successive SDP Relaxation Method (SSDPRM) can be written as 
follows. 

Algorithm 1.1. (SSILPRM) 

Step 0: Choose a D C D. Let k = 0. 

Step 1: If Ck = (the convex hull of F), then stop. 

Step 2: Let C k+1 =ft(C k ,D ). 

Step 3: Let k = k + 1, and go to Step 1. 

Algorithm 1.2. (SSDPRM) 



Steps 0, 1 and 3: The same as the Steps 0, 1 and 3 of Algorithm 1.1. 
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Step 2: Let C k+1 = N+(C k , D ). 

To connect these algorithms to the Lovasz-Schrijver procedures, we need to introduce 
some additional notation. For every pair of closed convex cones /C and T in R 1+n , define 



M{K,T) 
M + (K,T) 



v (X Xv T \ A>0 «€^, V€5», 

V ^ ^ v / /y«, > o, Vv e V, Vu> e K* 



X Xv T \ ^ - 0> v e C , V e S n , Y e S]+ n 
\ W : Vi = Vih ie {m + 1 ,m + ' 2 ,---,n}, 

' v T Yw > 0, Vv e V, \/w e K* 

Let To and /Co be closed convex cones given by 

Tl = c.cone {{ ( a{Co ^ d) I € : d G A, 



> . 



/Cn 



A 



a(C ,d) , ^_ r)1+ , ) 



Av 



G R 1+n :i)6Cn, A > 



(Note that To itself is defined as the dual of To-) If C C C is a compact convex relaxation 
of F and 

JC=U , A 1 e i2 1+m : » e C, A > ol , 



Xv 



then 



M{C,D Q ) 



veR n :(\ v ^ G A4(/C, T )j 
veRn '- { I v) zM+fcTo) 



Algorithms 1.1 and 1.2 specialized to (2) with Vf = { v i— v i, — v i+ v i, i £ {m + l,m + 
2, . . . ,n}} can be stated in the following forms, which are essentially the Lovasz-Schrijver 
procedures. 

Algorithm 1.1H (Homogeneous form of Algorithm 1.1) 
Step 0: Choose a D C D. Define To and /C as above. Let k = 0. 



Step 1: If /Cfc = c.cone 



v 



: v E F> ] then stop. 



Step 2: Let /C fc+ i = {Fe : Y G M(JC k , T )}. 
Step 3: Let k = k + 1, and go to Step 1. 



Algorithm 1.2H (Homogeneous form of Algorithm 1.2) 
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Steps 0, 1 and 3: The same as Steps 0, 1 and 3 of Algorithm 1.2H, respectively. 
Step 2: Let JC k+1 = {Ye : Y G M+()C k , To)}- 

In this paper ej denotes the jth unit vector and e denotes the vector of all ones (the 
dimensions of the vectors will be clear from the context). 

2 SCRMs applied to LCP with an a priori bound 

Let M e R txt , q e i?^ be given. Consider the LCP in the following form. 

(LCP) Find x, s such that Mcc + q = s, 

x > 0, s > 0, 

= 0, Vi e {1,2,...,£}. 

buppose we are given B(i,r) = {ue R 2£ : ||u - £|| < r}, an Euclidean ball containing 
a solution of the LCP. (In the case of rational data (M, q), we can take B centered at the 
origin with the radius bounded above by a polynomial function of the "bit size" of the data 
(M, q).) For the rest of this section, we assume that the Euclidean ball with center £ = 
and the radius r (r is assumed given) contains some solution of the LCP. 

Under the boundedness assumption above, it is particularly easy to model any LCP as 
a 0-1 mixed integer programming problem, since the only nonlinear constraints of LCP can 
be expressed as 

Xi = 0, or Si = 0, e {1,2,...,£}. 

Balas' method [1] can be directly applied to such formulations. We can also apply some 
variants of the Lovasz-Schrijver procedures [8] to the mixed integer programming feasibility 
problem: 

Find x, s and z such that M x + q = s, 

< x < rz, < s < r(e — z), 

ze{0,i} £ . 

Note that we can eliminate the variable vector s from the formulation and apply the 
SSILPR and SSDPR Methods to the following formulation: 

< Mx + q < r(e- z), 
< z < e, < x < rz, 

zf-Zi < 0, -zf + Zi < 0, ie{l,2,...,£}. 

To apply the SCRMs, we can take 

Co = \ V ~ [z ) E R • < z < e, < x < rz ) ' 

m = £, n = 2£, 

= {(vf-Vi), (-Vi+Vi), ie{m + l,m + 2,...,n}}. 
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Both algorithms, SSILPRM and SSDPRM presented in Section 1, terminate in at most 
I steps. This fact can be proved easily, using the results of Balas [1], Sherali and Adams 
[10], Lovasz and Schrijver [8], or Kojima and Tungel [5, 6]. For computational experience on 
similar algorithms for similar problems see [3], [12]. In the next section, we give the details 
of a proof of such a convergence result when the methods are applied to a formulation of 
Pardalos and Rosen [9]. 



3 SCRMs applied to Pardalos- Rosen formulation of 
LCP 



We will illustrate the convergence proof on a formulation of (LCP) by Pardalos and Rosen 
[9]. They homogenize the vector q with a new continuous variable a, then they maximize 
a. 

(MIP a ) maximize a 

subject to < Mx + qa < e — z, 

< x < z, < a < 1, z e {0, if. 

/ a 

Note that x \ = is feasible in (MIP a ) and, it is easy to see that (MIP a ) has an optimal 



solution with a* > iff the (LCP) has a solution (or solutions) [9]. Moreover, if 

x* 



(a*\ 

x* 
\z*J 



IS 



an optimal solution of (MIP a ) with a* > then — solves the (LCP) [9]. One advantage 

a* 

of (MIP a ) is that it does not require the introduction of large, data dependent constants 
(such as r in the previous section) or their a priori estimates. Now, we take 









c = < 


v = 


X 






w 



m = 



e R 



£ + 1, n = 2£ + l, 



1+21 



< Mx + qa < e - z, 
< x < z, < a < 1 



i € {m + 1, m + 2, 



n}} 



We have an analog of a very elementary but also a key lemma (Lemma 1.3 of [8]) of 
Lovasz and Schrijver (and their proof technique is adapted here). In what follows, we refer 
to the vectors in the space of fCk by v. At the same time, we refer to different subvectors 
of v by different names, such as x, a etc., to keep the correspondence of elements of v and 
the original formulation of F clearer. The proof of Lemma 1.3 of [8] leads to the following 
analogous result in our case. 

Lemma 3.1. Let D D {±e m+ i, ±e m+2 , . . . , ±e„}. Then the sequence of convex cones 
{K>k : k > 0} given by Algorithm 1.1H satisfies 

JCk+i C (JC k n{v.Xi = 0}) + (JC k n {v : (Mx + qa)i = 0}) , 

for every i G {1,2,..., £}, and for every k > 0. 
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Proof: Let w 



/1\ 

a 
x 



G KLk+i- Fix j G {1, 2, . . . ,£} arbitrarily. By the definition 



of -D and To, the unit vector e is in To- Hence, by the definition of .M(/Cfc,To), 
K-k+i ^ K-k for every A > 0. Therefore, w G /Cfc. If = or (.Ma; + qro;)j = then the 
statement of the lemma clearly holds. So, without loss of generality, we assume Xj > 
and (Mx + qa)j > 0. Let Y G Al(/Cfc,To) such that w = Ye . By our choice of the 
cone To, we conclude that Ye n+ j and Y(e — e n+ j) are both in /C&. Note that 



w 



w + w, 



where w 



Ye n+ j and w 



Y(e c 



-n+j 



). We will refer to the x and z parts of the vector 



w by x, z etc. (Similarly for w.) First, since by the definition of Ai{K,k,To), Vi = Vu 
for every i G {m + 1, m + 2, . . . , n}, we have 5j = which implies Xj = 0. Therefore, w 
lies in the cone (/Co H : Xj = 0}). Second, since Xj > 0, must be positive. Therefore, 
(l/zj)w G /Cq. Since ^ = Vjj for every % G {m + 1, m + 2, . . . , n}, i 



2j. 



So, 



with its .Zj entry equal to 1. Thus, (iVfa; + qfa)j = 0. Hence, w is in the cone 

(fCk H : (iVfa; + qr«)j = 0}). Since the argument above is independent of the index j 

the proof is complete. ■ 

Note that the conclusion of the above lemma also applies to the SSDPR Method since 
SSDPR Method yields at least as tight relaxations as the SSILPR Method. 



Theorem 3.2. Both algorithms, Algorithm 1.1H and 1.2H terminate in I iterations when 
applied to the formulation (MIP a ) with our choice ofVF, Co and D above. 



Proof: First note that 

c.hull(F) C <( 




/1\ 

a 

x 

w 



G/G 



> , VA; > 0. 



Next, let i, j G {1, 2, ...,£}, i ^ j. Since x > and Mx + qa > 0, for all v G /C^, for 
every A > 0, 

[(/Cjb n {w : a;i = 0}) + (/C fc n {v : (Ma; + qra)^ = 0})] n {v : Xj = 0} 
= (/C fc n{v.Xi = 0, Xj = 0}) + (/C fe n : ^ = 0, (Mx + qa)i = 0}) . 

Similarly, for the intersection with {v : (Mx + qa)j = 0} . Now, we apply Lemma 3.1 
repeatedly to conclude that /Q is the homogenization of the convex hull of all solutions 
of the LCP that lie in the original relaxation Cq. i 
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4 SCRMs applied to the smaller formulation of LCP 
with explicit treatment of disjunctive constraints 



Now, we consider a formulation with fewer variables and constraints. 



(LCP a ) maximize a 

subject to Mx + qa > 0, x > 0, a > 0, 
e T (M + I)x + {e T q + l)a < 1, 
Xi(Mx + qa)i = 0, i G {1, 2, . . . , £}. 



It is easy to see that 



is feasible in (LCP a ), and it is also easy to observe that 



(LCP a ) has an optimal solution with a* > iff the (LCP) has a solution(or solutions). 

fx*\ x* 
Moreover, if ( x \ is an optimal solution of (LCP a ) with a* > then — solves the (LCP). 

Note that the inclusion in Lemma 3.1 can sometimes be strict for the SSILPR and SSDPR 
Methods. 

We explicitly include the variable vector s in our discussion in this section, for the sake 
of presentation. Let 







fx) 




c = < 


v = 


s 


e R 2£+1 




< 


{a) 





s = Mx + qa > 0, x > 0, a > 0, 
e T (M + I)x + (e T q + l)a < 1, 



In this section, we will describe another Successive Convex Relaxation Method based on the 
ideas of Balas [1], Lovasz and Schrijver [8]. This method will use only Linear Programming 
(LP) relaxations. We describe the method in the original space of F and C . Let T(C ) 
denote the set of facet defining inequalities for C . J-(C ) is the input of the algorithm 
which we introduce now. 



Algorithm 4.1. Step 0. k = 0. 
Step 1. T{C k+l )=T{C k ). 
Step 2. For every inequality 

i 

- ^2 {uiXi + u l+i Si) - u 2 i + ia < u 

i=l 

in T(Ck) and every j € {1, 2, . . . , £} solve the LP problems 

(Pj) minimize 

subject to = 1, &j = 0, $ U) e )C k , 

and 

{Pi+j) minimize u T ^ + ^ 

subject to = dS J) = 1- t (e+j) < K k . 
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If (Pj) is infeasible then add the equation Xj = (or the inequality Xj < 0, since the 
inequality Xj > is already included) to T{Ck+\)- If (Pt+j) is infeasible then add the 
equation sj = to F(Ck+i)- Otherwise, let (£^)* and denote the optimal solutions 

of {Pj) and (P l+j ) respectively Define ^ = Uj - w T (£ (i) )*, y e+j = u l+j - u T (£ ii+j) )* . Add 
the inequality 

- ( U i X i + U l+i S i) ~ Vj X j ~ Vl+j S j - U 2 £+ia < U 

to T(C k+ i). 

Step 3. Let k = k + 1, and go to Step 1. 

Note that in iteration k, the algorithm solves (2£\J 7 (Ck)\) LP problems. 

Theorem 4.2. Let C^, k e {1,2,...} be the sequence of convex relaxations generated by 
Algorithm J,..l. Then Ct = c.hull(F). 

Proof: We think of /C& for all k > 0, as a subset of R l+<y2l+l \ with the 0th compo- 
nent being the homogenizing variable, the next £ components representing x, the next I 
components representing s and the last component representing a. Note that 

/Ci C (/C n {v : Xj = 0}) + (/C n {v : s,- = 0}) 

iff 

1C\ D (/C* + {-e,-}) n (/C* + {-e, +i }) . (3) 

(We used the fact that /C C i?^_ + ^ 2f+1 \) Therefore, if we ensure the latter inclusion, then 
Theorem 3.2 applies and we can conclude the convergence of the method in I iterations. 
Recall that every vector u e /Cq represents a valid inequality 

1 

- ^2 {UiXi + u l+i Si) - u 2 i + ia < u 

2 = 1 

for Co. To ensure the inclusion (3), it suffices to prove: 

"For every u, w G /Cq such that Ui = Wi, Vz ^ {j, £ + j}',Uj > Wj, ui + j < wi + j, 

we have y e /C*, where Vi = Ui,Wi ^ j; yj = Wj." 
This is equivalent to proving the fact that if the two inequalities 

£ 

- ( u i x i + u l+i Si) - u 2e+1 a < u , and 

2 = 1 

- ^ (UiXi + - lOj-Xj - w i+j Sj - u 2 i + ia < u 
are valid for Co, then 

- ( U i X i + U l+i S i) ~ W i X i ~ U t+jSj - U 2t+1 a < U 
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is valid for C\. To compute all such inequalities defining Ci, we solve for every valid 
inequality 

£ 

- ^ (UiXi + U l+i Si) - U 2 l+lOi < Uq 
i=l 

for C and every j € {1, 2, . . . , £}, the linear programming problems 

maximize 3 

subject to @ej + Sei + j u, 

and 

maximize 7 

subject to nej + jei + j u. 

Here, denotes the partial order induced by the convex cone /Q (that is, u 1 u 2 
iff {u 2 — u 1 ) e KLq). Note that both problems are always feasible. Therefore, each of them 
either has an optimal solution or is unbounded. If both LPs have optimal solutions, say 
j5* and 7* then we set Wj = Uj — j5* and ui+j = ui + j — 7*. Since the above two problems 
are LPs, we can equivalently solve their duals. Namely, we solve the LPs: 

(Pj) minimize u T £^ 

subject to = 1, &j = 0, ^ (J) e /Co, 

and 

(Pi+j) minimize u T ^ l+ ^ 

subjcet to Cf +J) = 0, dS j) = 1, G /C . 

These latter two linear programming problems are precisely the ones used by Algorithm 
4.1. Notice that since their duals are either unbounded or have optimal solutions, these 
LP problems either have optimal solutions or are infeasible. When (Pj) is infeasible, the 
equality Xj = is valid for F and the algorithm adds this equality to the describing 
inequalities of C&. Similarly, when (Pi+j) is infeasible, Sj = is valid for F and the 
algorithm behaves correctly in this instance. (In either instance, the inclusion (3) is 
obviously satisfied for j.) However, the proof is not yet complete; because, the arguments 
so far ensure the inclusion (3) when the algorithm is ran for every valid inequality of Co. 
So, next we prove that what the algorithm does (using only the facets of C ) suffices. To 
see this, we need to prove that to derive the facets of JCi, it suffices to start with a facet 
u of /C in the above procedure. Suppose u, w e /Q satisfy the above conditions but u is 
not facet inducing for KT Q . (We will prove that the valid inequality derived from u and w 
is implied by some other inequalities derived from some facets u 1 , u 2 , . . . , u l of /Co.) Since 
u is not facet inducing for /C , u is not an extreme ray of K,* Q . Hence, there exist extreme 
rays w 1 , w 2 , . . . ,u e of /Cq such that for some A r > 0, r € {1, 2, . . . , £}, X)f=i X = 1 the 
following conditions are satisfied: 

£ 

u = ^2 Ku r , u r = u , Vr G {1,2,..., £}. 

r=l 
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Note that u r is facet inducing for each r. Let £ r be the optimal solution of (Pj) above for 
the objective function vector u r . Let £* be an optimal solution of (Pj) when the objective 
function vector is u. We claim that there exists £ € /C such that 

K) T £ = K) T r, Vr€{l,2,. ..,*}, £ = 1, & +j = o, leK,. 

(This claim follows from Farkas' Lemma, using the facts that u r G /Q, Vr and £ r 6 /C , Vr.) 
Thus, we have 

EAr(« r ) r r = ^ > « T r- 

r=l 

Therefore, the inequality obtained from u is equivalent to or dominated by a nonnegative 
combination of the inequalities obtained from u r which induce facets of /Co. The proof is 
complete. | 

We illustrated a derivation and convergence proof for a successive relaxation method 
(closely related to Balas' approach and analogous to a suggestion of Lovasz and Schrijver 
[8]) based on Lemma 3.1 and Theorem 3.2. Algorithm 4.1 is an analog of a method based on 
relaxations Nq(JC) from [8] (which is concerned with the case of 0-1 integer programming). 
For the relationship of the methods of [1] and [8], see Balas, Ceria and Cornuejols [2]. (Balas' 
method [1], in essence, corresponds to defining 

JC k+1 = (JC k n {v : x k+l = 0}) + (/C n {v : (Mx + qa) k+l = 0}) .) 

Let C^ k A \ k > denote the projection of C k generated by Algorithm 4.1 onto the coor- 
dinates {^j- Let k > denote the projection of C k , generated by Algorithm 1.1, as 

used in Section 3, onto the coordinates ( X | . Let /Cl 4 ^ denote the convex cone associated 



with C^K From the proof of Theorem 4.2, it is easy to see that 

4?i = n [(4 4) n {v : Xi = 0}) + (4 4) n{v:s z = 0})] . 

i=l 

Therefore, the proofs of Theorems 3.2 and 4.2 imply that 

if C ( o ] D C^ 3) then C ( k 4) D C [ k 3) for all k > 0. 

Thus, the SSILPR Method (Algorithm 1.1) as applied in Section 3 to (MIP a ) converges at 
least as fast as Algorithm 4.1 applied to (LCP Q ). 



5 SCRMs applied to the smaller formulation of LCP 
with an implicit treatment of the disjunctive con- 
straints 

We have already seen various ways of applying SCRMs to LCP problems. Since the methods 
proposed in [5, 6] only require a formulation of the feasible solutions by quadratic inequali- 
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ties, we are also interested in applying the methods of [5, 6] to the following formulation: 

fa\ t+l Mx + qa > 0, x > 0, a > 0, 1 
UJ : e T (M + I)x + (e T q + l)a < 1 j' 

and 

V F = {xi(Mx + qa)i < 0, i E {1, 2, ...,£}} . 

The general theory of Kojima-Tuncel [5] implies that their SSDPR and SSILPR Methods 
converge. It would be interesting to characterize the conditions under which the Algorithms 
3.1 and 3.2 of [6] converge in at most I iterations for the above description of Vf and C . 
Also see [7], where the authors derived some necessary and some sufficient conditions for 
the finite convergence of SCRMs. 




6 A general linear complementarity problem 

Let A : R £ — > R l , a linear transformation, q e R l and K, C R l a pointed, closed convex 
cone with nonempty interior, be given. Consider the Complementarity Problem (CP): 

(CP) Find x, s such that A(x) + q = s, 

x e /C, s e /C*, (x,s) = 0, 

where /C* is the dual of /C: 

/C* = {seR 1 -. (x,s) > 0, VxelC}. 

Since /C is a pointed, closed convex cone with nonempty interior, so is JC*. Such problems 
were studied recently, in the context of interior-point methods [11]. We pick r\ € int(/C), 
f] G int(/C*) and we can solve instead the optimization problem 

a 

x £ JC, [A(x) + qa] e JC*, a > 0, 
{fj, x) + (r), A{x) + qa) + a < 1, 
(x,A(x) + aq) = 0. 



x e JC, [A(x) + qa] e JC*, a > 0, 1 
(ff, x) + (rf, A(x) + qa) + a < 1 J' 

Note that Co is always a compact convex set (see the next theorem). We also pick 

V F = { {x,A(x) + aq), -(x,A(x) + aq) j. 

Theorem 6.1. (i) C is a compact convex set. 
(ii) (CP a ) has an optimal solution with a* > iff (CP) has a solution (or solutions). 



(CP a ) maximize 
subject to 



We choose 



Co = <[ a ) e R i+1 

X 



12 



(Hi) If y^,* ) i> s an optimal solution of (CP a ) with a* > then the pair of vectors 
'x* 1 



(^,—A(x*) + q) solves (CP). 
\a* a* J 



Proof: 

(i) We only need to show that C is bounded; because, C is a closed and convex subset of 
R e+1 by definition. Assume on the contrary that we can take an unbounded direction 

(£)#Oi„C„; 

( Ax ) ^ °' Ax e /C ' Aa - °' ^ Ax ^ + qAa ^ e 
(f),Ax) + (r),A(Ax) + qAa) + Aa < 0. 

Since each term in the left hand side of the last inequality is nonnegative, we have 

(fj, Ax) = and Aa = 0. 



Since rf € int(/C*) and Aa; e /C, the first identity above implies that Aa; = 0. Thus, 

Aa Wo. 

J 



we have a contradiction to 



Aa; 



(ii) Suppose (CP Q ) has an optimal solution with a* > 0. Then x = — e /C, 

V a; / n 1 

s = — A(x*) + q E IC*. We have 

a* 

(a;,s) = (a;,^l(a;) + q) = j^(x*, A(x*) + a*q) = 0. 

Therefore, (a;, s) solves (CP). For the converse, let (x,s) be a solution of (CP). Let 

1 x 
( = {fj, x) + (rf, s) > 0, a* = and x* 



C + i C + i' 

Then ^* ^ i s a feasible solution of (CP Q ). But the feasible region of (CP Q ) is compact 

and nonempty, its objective function is linear, hence, (CP Q ) has optimal solution (or 
solutions). Since we already showed a solution with positive objective value, the 
optimum value is positive. 

(iii) This claim follows from the proof of (ii). 



Theorem 6.1 shows that we can apply SCRMs to (CP Q ) with the above C and Vf and 
solve the original, general problem (CP). 



13 



References 



[1] Balas, E. (1974), Disjunctive programming: Properties of the convex hull of feasible 
points. Management Science Research Report 348 GSIA, Carnegie Mellon University, 
Pittsburgh, PA, USA. 

[2] Balas, E., S. Ceria and G. Cornuejols (1993), A lift-and-project cutting plane algorithm 
for mixed 0-1 programs. Mathematical Programming 58 295-323. 

[3] Ceria, S. and G. Pataki (1998), Solving integer and disjunctive programs by lift-and- 
project, Proceedings of the Sixth IPCO Conference. 

[4] Cottle, R. W., J.-S. Pang and R. E. Stone (1992), Linear Complementarity Problem, 
Academic Press. 

[5] Kojima, M. and L. Tungel (1998), Cones of matrices and successive convex relaxations 
of nonconvex sets. SIAM J. Optimization, to appear. 

[6] Kojima, M. and L. Tuncel (1998), Discretization and localization in successive convex 
relaxation methods for nonconvex quadratic optimization problems. Technical Report 
B-341, Dept. of Mathematical and Computing Sciences, Tokyo Institute of Technol- 
ogy, Tokyo, Japan. Also: Research Report 98-34, Department of Combinatorics and 
Optimization, University of Waterloo, Waterloo, Ontario, Canada. 

[7] Kojima, M. and L. Tungel (1999), On the finite convergence of successive SDP re- 
laxation methods. Technical Report B-354, Dept. of Mathematical and Computing 
Sciences, Tokyo Institute of Technology, Tokyo, Japan. Also: Research Report 98-36, 
Department of Combinatorics and Optimization, University of Waterloo, Waterloo, 
Ontario, Canada. 

[8] Lovasz, L. and A. Schrijver (1991), Cones of matrices and set-functions and 0-1 opti- 
mization. SIAM J. Optimization 1 166-190. 

[9] Pardalos, P. M. and J. B. Rosen (1988), Global optimization approach to the linear 
complementarity problem. SIAM J. Scientific and Statistical Computing 9 341-353. 

[10] Sherali, H. D. and W. P. Adams (1990), A hierarchy of relaxations between the con- 
tinuous and convex hull representations for zero-one programming problems, SIAM J. 
Discrete Mathematics 3 411-430. 

[11] Shida, M., S. Shindoh and M. Kojima (1997), Centers of monotone generalized com- 
plementarity problems. Mathematics of Operations Research 22 969-976. 

[12] Stubbs, R. A. and S. Mehrotra (1996), A branch-and-cut method for 0-1 mixed convex 
programming, Mathematical Programming, to appear. 



14 



